Short-ranged and short-lived charge-density-wave order and pseudogap features in 

underdoped cuprates superconductors 
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The pseudogap phase of high-Tc cuprates is controversially attributed to preformed pairs or to 
a phase which coexists and competes with superconductivity. One of the challenges is to develop 
theoretical and experimental studies in order to distinguish between both proposals. Very recently, 
researchers at Stanford have reported [M. Hashimoto et al, Nat. Phys. 6, 414 (2010); R.-H. He et 
al., Science 331, 1579 (2011)] angle-resolved photoemission spectroscopy experiments on Pb-Bi2201 
supporting the point of view that the pseudogap is distinct from superconductivity and associated 
to a spacial symmetry breaking without long-range order. In this paper we show that many features 
reported by these experiments can be described in the framework of the t-J model considering 
self-energy effects in the proximity to a d charge-density-wave instability. 

PACS numbers; 74.72.-h, 74.25.Jb, 71.10.Fd, 79.60.-i 
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Although the solid state physics community agrees on 
the existence of the pseudogap (PG) phase in under- 
doped high-Tc superconductors, the origin of the PG 
and its relation with superconductivity (SC) is not clear 
from experimental and theoretical studies. Interpreta- 
tions run from descriptions where the PG is intimately 
related to SC to others where the PG is distinct from 
SC and both phases competcji Thus, experimental and 
theoretical studies are of central interest for solving this 
puzzle. 

Angle-resolved photoemission spectroscopy (ARPES) 
is a valuable tool for studying the PG phascj^ In the 
pseudogap phase ARPES experiments show, in the nor- 
mal state below a characteristic temperature T*, Fermi 
arcs^"- (FAs) (centered along the zone diagonal) instead 
of the expected Fermi surface (FS). Despite the general 
consensus about the existence of FAs, their main char- 
acteristics are controversial. Some experiments^ii^ report 
that near the antinode exists a single gap which is nearly 
independent of temperature. Moreover, in the supercon- 
ducting state the gap follows the expected d-wave behav- 
ior along the FS. In contrast, other experiments^!^ show 
that the gap follows the d-wave behavior near the node 
but deviates upward approaching the antinode. In our 
opinion the difference in the spectra reported by differ- 
ent groups can not only be attributed to the different 
systems under study. For instance, similar samples at 
similar doping (see Refs.^] and ^]) also show the men- 
tioned discrepancy. Although the reason for these differ- 
ences is not clear, results of Refs.fS] and 3 suggest that 
FAs are tied to thepreformed pair scenario, while results 
from Refs. Q and @ claim that arcs are associated to an 
order which is distinct from SC. 

Recently, ARPES results on Pb-Bi2201 were reported 
(Refs.Q and ^]) showing that the PG phase is distinct 
from SC and associated to a spacial symmetry breaking 
without long-range order. Since these experiments chal- 
lenge scenarios about the origin of the PG and FAs, their 
theoretical support is of high and broad interest. 
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FIG. 1. (Color online) Energy distribution curves in the nor- 
mal state from a high temperature (left panel) to a lower 
temperature (right panel) near the d-CDW instability along 
the antinodal cut shown by an arrow in the inset of the left 
panel. 



In the framework of the t-J model at mean-field level 
the carriers show a band dispersion efe = — 2ie//[cos(fc2:)-l- 
cos(fcy)]— 4tg^j cos(fci,) cos(A:j,)— /I, with tg// = (t| + rj) 



and t' 



eff 



i'§. i, t' and J are the hopping between 
nearest-neighbor, next-nearest-neighbor and the nearest- 
neighbor Heisenberg coupling, respectively. S is the dop- 
ing away from half-filling, fi is the chemical potential, 
and rJ is the J-drivcn hopping term. This mean-field 
homogeneous Fermi liquid (HFL) phase presents an in- 
stability to a d charge-density-wave (d-CDW) stated or 
a flux phase^'^ (FP) below a characteristic temperature 
Tp p which decreases with increasing doping (see Ref. [TTj 
and references their in). 

Beyond mean-field level, in the HFL phase and in 
the proximity to the d-CDW instability, the self-energy 
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FIG. 2. (Color online) Intensity maximum dispersion, for sev- 
eral temperatures in the normal state, taken along the antin- 
odal cut. Solid line is the predicted mean-field dispersion. At 
mean-field level the effective parameters are teff/t ~ 0.08 
and t'^ff/t ^ —0.017, while the effective parameters that 
reproduce the effective band at high temperature (dotted 

line) are t^ff/t ~ 0.21 and t'^ff/t ^ 0.03. Note that the 

Fermi wave vector kp is different from the mean-field value 
k^^ ~ (7r,0.l7r). 

S(k,aj) readsii 

/mE(k,w) = --^^7^(q,k)/TOx(q,w - tk~q) 

' q 

X [nF(-e/e-g) + rts(t^ - e/c-g)] , 

where uq and np are the Bose and Fermi factors, re- 
spectively, and A^s is the number of sites. xCQj'*^) = 
(f)2[(8/J)r2 -n(q,w)]-i is the d-CDW susceptibility 
that diverges, at w = and q = (7i',7r), and at Tpp. 
n(q, lS) is the electronic polarizability calculated with a 
form factor 7(q, k) = 2r[sin(fca; — g^/S) — sin(fcj^ — qj,/2)]. 
Note that since q ^ (tt, tt) the form factor 7(q, k) trans- 
forms into ^ \cos{}ix) — cos(fcj,)] which indicates the 
d-wave character of the d-CDW. Using the self-energy 
S(k, w), the spectral function is calculated as usual. This 
scenario leads to a PG and FAs which depend on dop- 
ing and temperature in close agreement with the experi- 
ment, and occur via the coupling between quasiparticles 
(QPs) and the c?-CDW soft mode in the proximity to the 
instabilityJii^i^ Here after t' jt = -0.35 and J jt = 0.3. 
The hopping t and the lattice constant a of the square 
lattice are used as energy and length units, respectively. 

We show here that many aspects observed in Re fs.M 
and @ can be discussed in the framework of Refs.[ll|. 
[l^ . and Ts'] showing that a theory where the PG phase 
is distinct from SC and associated with short-range fluc- 
tuations in the proximity to a symmetry breaking insta- 
bility is consistent with the experiment. 

Similar to Fig. 1 g-1 of Ref. 0, in Fig. 1 we present 
energy distribution curves (EDC) for momenta along the 
antinodal cut indicated in the inset. With no loss of gen- 
erality, results are given in the HFL phase for doping 
6 = 0.10 and for several temperatures from T/t = 0.035 
to T/t = 0.019 close but above Tpp/t - 0.018. At high 
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FIG. 3. (Color online) (a) Energy distribution curves for mo- 
mentum (tt, 0) for several temperatures, (b) The same as (a) 
for momentum kp. The inset shows the energy position of 
the maximum of the energy distribution curves as a function 
of temperature for k = (7r,0) (circles) and kp (triangles). 



temperature (T/t = 0.035), EDC resemble a parabolic- 
like dispersion as expected for usual metals. Lowering 
temperature, a PG like feature becomes evident near the 
antinode and, at the same time, a back-bending occurs in 
the dispersion. For summarizing the main results_of Fig. 
1, in Fig. 2, we show (similar to Fig. In in Ref.|7j]) the 
intensity maximum dispersion for several temperatures. 
At high temperature {T/t ^ 0.1), no PG is observed near 
the antinode and a parabola-like dispersion (dotted line) , 
which crosses uj — at hp ^ (tt, O.IGtt), is obtained. In 
the figure, we have also plotted results for the mean-field 
band (solid line). With decreasing temperature a PG 
opens near the antinode and, as in the experiment, a 
back-bending (marked by an arrow) is observed beyond 
kp, in addition, the energy of the band bottom {Epot) at 
(tt, 0) moves deeply to lower binding energies. In Ref.0, 
the existence of a back-bending beyond kp was discussed 
as a result that rules out preformed pairs and favors a 
scenario where the PG is distinct from SC. Present re- 
sults are consistent with that interpretation as discussed 
below. 

Following the comparison with results in Ref.[3], in 
Figs. 3(a) and 3(b) we have plotted EDC at (7r,0) and 
kp, respectively, for several temperatures. For k = (tt, 0), 
the energy of the maximum of the EDC moves, with in- 
creasing temperature, toward Epot as shown by circles in 
the inset. Figure 3(b) shows that, with increasing tem- 
perature, the PG at kp closes and, at the same time, 
fills as in the experiment. In the inset the energy of the 
maximum of the EDC at kp is plotted (triangles) as a 
function of temperature showing that, for T/t > 0.06 no 
PG is observed. Since EDC at kp show a maximum at 
cj = for T/t > 0.06, we identify the pseudogap tem- 
perature with T* /t ~ 0.06. In addition, we also identify 
the dotted line in Fig. 2 with the effective metallic band 
(recovered at high temperature) discussed in Ref.Q- Fi- 
nally, we note that in the present theory the opening of 
the PG is not associated to an abrupt transition but to 
a smooth crossover. 

Using the accepted value t = 0.4 eV, T* - 270 K, 
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FIG. 4. Intensity maximum dispersion for several cuts, par- 
allel to the antinodal cut, at T/t = 0.020. Dashed line in 
each panel is the predicted mean-field dispersion for the cor- 
responding cut. 



Esot 25 meV, and the energy at the back-bending ^ 40 
meV. Since our description does not require phenomeno- 
logical parameters as a value for the PG energy and its 
temperature dependence, as in other approaches'^ the 
agreement with the experiment can be considered satis- 
factory. 

Next, we show predictions of our theory for other cuts 
on the Brillouin zone. Similar to Fig. 2, in Fig. 4 we have 
plotted, for T/t — 0.020, the intensity maximum disper- 
sions for cuts parallel to the antinodal cut. In each cut, 
we fix to be tt, O.Ott, O.Stt, O.Gtt, and O.Stt, and vary ky. 
While for cuts close to the antinode [(tt, fcy), (O.Qtt, fcj^), 
and (O.Stt, ky)] a PG and a back-bending are clearly ob- 
served, for cuts near the node [(O.Gtt, fcy) and (O.Stt, fcj^)] 
the back-bending and the PG wash out and a FS crossing 
occurs. This behavior is consistent with the prediction of 
FAs by present theory.—"— In addition, different to the 
cuts near the antinode (see Figs. 1 and 3), for cuts near 
the node the QP peaks (not shown) become sharper and 
the renormalized dispersion is closer to the bare mean- 
field band (dashed line) . This fact reflects that the effect 
of the d-CDW mode diminishes approaching the node.^^ 

It is worth to mention that a generic long-range CDW 
phase may show, qualitatively, some of the features pre- 
sented here such as, for instance, the back-bending be- 
yond kp. Some differences with the long-range CDW 
point of view are: (a) While in a CDW phase the transla- 
tional symmetry is broken, in the present theory, the sys- 
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FIG. 5. (Color online) (a) Energy distribution curves along 
the antinodal cut for T/t = 0.1. Thick solid blue line is the 
energy distribution curve for k = (tt, 0). Dashed red lines rep- 
resent a guide for the eye showing that the double peak struc- 
ture at both sides of a; = near k — (tt, 0) merges into only 
one peak away from k — (7r,0), which follows the predicted 
mean-field result. Note that at kp (dashed-dotted green line), 
no pseudogap is observed for this temperature, (b) Intensity 
maximum dispersion (solid line), the effective band (dashed- 
dotted line) discussed in text, and the mean-field dispersion 
(dashed line) along the antinodal cut for T/t — 0.1 (c) The 
same as (b) for the cut (O.Stt, fcy). 

tem is always in the HFL phase but, under the influence 
of short-range fluctuations, (b) While sharp QP peaks 
are expected in a CDW phase, our results show broad 
structures as in the experiment. These broad structures 
are intrinsic and not related to the instrumental resolu- 
tion as suggested in Ref.[14]. For explaining the broad 
spectra observed in the experiment, in Ref.Q a static 
density wave order with a finite correlation length was 
assumed. In addition, it was suggested that the observed 
fast sink, with decreasing temperature, of the band bot- 
tom implies an incommensurate character of the charge 
order. Both, the broad spectra and the temperature be- 
havior of the band bottom are naturally described by the 
present picture. We think that this fact suggests, besides 
of a short correlation length, a slowly dynamical charac- 
ter of the charge order, and these conditions are fulfilled 
in the fluctuating regime near the instability. Since the 
discussion about the interpretation of the PG as a phase 
with long-range order versus a phase without long-range 
order is highly controversial (see Ref.jlHl for a review), we 
expect that our finding may contribute to the distinction 
between both scenarios. 

Now, we discuss some characteristics of the high tem- 
perature phase. In Figs. 5(a) and 5(b), we present EDC 
and the corresponding intensity maximum dispersions. 
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FIG. 6. (Color online) EDC obtained for various supercon- 
ducting gaps values A (see text). Inset: Energy position of 
the leading gap of the EDC as a function of the FS angle (j) 
for the normal state (dashed line) and in the superconducting 
state (solid line). 

respectively, for the antinodal cut at T/t = 0.1. Note 
the broad feature with structure at both sides of a; = 
at (tTjO) [thick solid blue line in Fig. 5(a)]. This double 
peak structure merges into only one feature away from 
(tt, 0) approaching the predicted mean-field band [dashed 
line in Fig. 5(b)]. The parabola-like band, discussed in 
Fig. 2 at high temperature, follows the effective band 
[dashed-dotted line in Fig. 5(b)] for a; < 0. We note 
that the closing of the PG at T/t ^ 0.06 is mainly due 
to a filling effect which occurs from the merging of the 
broad peak at w < with the one at w > 0. In Fig. 
5(c), we have plotted, for the cut (O.Stt, ky), the intensity 
maximum dispersion of the EDC at T/t — 0.1. Different 
to the antinodal cut, the intensity maximum dispersion 
is closer to the mean-field band (dashed line) for both, 
w < and w > 0. From the experimental point of view, 
this result means that above T*, a unique set of param- 
eters might not describe the dispersion for all cuts, as it 
should be if the high temperature phase can be associ- 
ated with a usual metal. It is interesting to see whether 
this prediction may be checked by the experiment. 

In Ref.[7|, it was reported that with decreasing tem- 
perature, the spectra and the leading edge gap near 
the antinode smoothly become broad and move toward 
higher binding energy, respectively, with no significant 



change below the superconducting critical temperature 
Tc- At first sight, it would expect that d-CDW fluc- 
tuations would be washed out once superconductivity 
sets in. However, considering that the superconduct- 
ing gap is smaller than the pseudogap, we can approxi- 
mate the Green function in the superconducting state as 
Gg^{k,iuj) = G~'^{k,iuj) + AlG{k, -iuj), where G{k,iuj) 
is the Green function describing the spectral functions of 
Fig. 3(b), and — A(cosfca; — cosfcy)/2 is a d-wave 
superconducting gap. In Fig. 6, we present EDC at kp. 
In the normal state (A = 0) (assumed to be the previous 
result for T/t — 0.025), the peaks at both sides oi u — 
are broad and asymmetric. With increasing A (decreas- 
ing temperature below Tc), the broad peak slowly moves 
to lower energies as in the experiment.^ It is interesting 
to note the presence of an inner shoulder which is a direct 
indication of two distinct energy scales below Tc. A sim- 
ilar shoulder was discussed in Ref. [7j, and more clearly 
reported in Ref. [8,] after improving the experimental res- 
olution. We think that an unambiguous experimental 
confirmation of an inner shoulder (or a peak) below Tc is 
of fundamental interest for solving the puzzle about one 
versus two energy scales. The reason for the absence of 
this shoulder in some ARPES experiments^!^ is an open 
question and out of the scope of the present paper. In 
the inset of Fig. 6 we plot the energy position of the 
leading gap of the EDC as a function of the FS angle (p 
from the antinode {(f) = 0) to the node (</> = 45 degrees). 
For A = (dashed line), the region near the antinode 
is gaped leading to Fermi arcs. For A/t — 0.07 (solid 
line) , the curve follows the d-wave behavior (dotted line) 
only near the node, and deviates upward approaching the 
antinode. It is interesting to see that for Pb-Bi2201, the 
gap follows this trend in similar way to Refs.[ij and [H 
(see Fig. 4S in Ref.iSj]). 

Concluding, we have shown that the scattering be- 
tween carriers and d-CDW fluctuations, which occur in 
the proximity to the d-CDW instability, describes many 
features observed in recent ARPES experiments. The 
proposed scenario supports the point of view that the 
pseudogap is distinct from superconductivity and associ- 
ated with short-range and short-lived CDW fluctuations. 

The authors thank H. Parent, R.-H. He, and R. Zeyher 
for valuable discussions. 
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